Solar Flare Intermittency and the Earth's Temperature Anomalies 
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We argue that earth's short-term temperature anomaUes and the solar flare intermittency are 
linked. The analysis is based upon the study of the scaling of both the spreading and the entropy of 
the diffusion generated by the fluctuations of the temperature time series. The joint use of these two 
methods evidences the presence of a Levy component in the temporal persistence of the temperature 
data sets that corresponds to the one that would be induced by the solar flare intermittency. The 
mean monthly temperature datasets cover the period from 1856 to 2002. 
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The historical recognition that the sun warms the earth 
has suggested a direct connection between the average 
global temperature and solar activity. Consequently any 
significant changes in solar activity should result in equiv- 
alent changes in the earth's global temperature. The lit- 
erature on the solar influence on the earth's temperature 
is quite extensive Q, |^ |^ , indicating the importance of 
the problem and that there are many issues that require 
further investigation. Herein we address the relation be- 
tween the statistics of solar flare activity and the fluctu- 
ations in the earth's global temperature. 

The dynamics of the sun's surface is turbulent, as is 
evidenced by changes in solar flare activity, with 11 or 
22 year solar cycles 0, 1^ and strong erratic fluctuations 
associated with solar flare intermittency ^ IE 111 - So- 
lar irradiance changes in accordance with the frequency 
of solar flares because each flare releases more energy 
than the background irradiance This time variation 
in the frequency of solar flares induces a similar pseudo- 
periodic cycle in the earth's average temperature, as well 
as produces trends that move the global temperature up 
or down for tens or even hundreds of years 0, 0, U- 
However, it is less evident that short-term (weekly and 
monthly) changes in the global temperature are tied to 
solar activity whose short-time fluctuations would have 
the intermittent dynamics of solar flares. Traditional 
measures, such as cross-correlation functions, would not 
show the connection between short-time fluctuations in 
the global temperature and solar flare activity, because 
of the strong nonlinear hydrodynamic interactions across 
the earth's surface. For example, the hydrodynamic in- 
teraction of the atmosphere over land and water, would 
suppress any direct correlation between the intermittent 
sun's irradiance and the earth regions' short-time re- 
sponse. 

This letter focuses on an alternate approach to estab- 
lishing the connection between the sun's irradiance and 
the earth's temperature fluctuations. A link between the 
two phenomena is detected through a detailed scaling 
analysis of the time series for the earth's temperature 
and the time series for the solar flare frequency. 

Considering a solar flare as an event, the time series for 



the number of solar flares has been interpreted as a wait- 
ing time distribution function between events. The so- 
lar flare waiting time distribution function is determined 
to be an inverse power-law probability density function 
(iplpdf ) 0, IE 01 • Such dynamical stochastic processes 
can be described by generalizations of random walks. A 
Levy-flight, for example, is such a process with a diverg- 
ing second moment. A Levy-walk, on the other hand, vis- 
its the same spatial sites as does a Levy-flight, but each 
step takes a finite time and the second moment is finite. 
The time necessary to complete a step in a Levy-walk is 
specified by an inverse power-law waiting time distribu- 
tion function, as first noticed by Shlesinger et al. It 
has been determined 0| that a Levy-walk can describe 
the intermittent solar fiare signal. Scafetta et al. 0, 0| 
established that, in general, the presence of a Levy-walk 
process in a given time series can be detected by the joint 
use of two separate scaling techniques, the Diffusion En- 
tropy Analysis (DEA) and Standard Deviation Analysis 
(SDA). We apply the same approach to the analysis of 
temperature data sets and compare them with the Levy- 
walk statistics induced by solar flare intermittency. 

We study the earth's temperature anomalies both glob- 
ally and locally, that is, North-South hemispheres and 
Land-Ocean regions, see Fig. 1. The technical term tem- 
perature anomalies denotes temperature departures from 
the 1961-1990 mean temperature value. These data are 
recognized in the geophysical community as among the 
most accurate data files for global air temperature, and 
global sea surface temperature (SST) |0. The basic data 
set of global earth's temperature anomalies (HadCRUT) 
is a combination of Land air temperature anomalies ^l| 
(CRUTEMl) and sea-surface temperature anomalies [IJ 
on a 5° X 5° grid-box basis. The merging of the two data 
sets is discussed in Parker et al. 12] and more recently 
in Jones et al. The Land surface time series are 

calculated using data from the Global Historical Clima- 
tology Network (Version 2) and sea-surface temperature 
anomalies from the United Kingdom MOHSST data set 
and the NCEP Optimum Interpolated SSTs (Version2). 
The data that we analyze are mean monthly tempera- 
ture anomalies downloaded from Climatic Research Unit, 
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FIG. 1: Global and local Earth temperature anomalies in 
Celsius degree (years: 1856-2002). 
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FIG. 2: Fig. 2a shows a computer-generated intermittent se- 
quences of impulses. Fig. 2b shows the frequency of impulses 
for 1700 time units. 



UK, [Ijl (Global and North-South temperatures, 1856- 
2002) and from National Climatic Data Center, USA, 
[15] (Land-Ocean temperatures, 1880-2002). 

The analysis of the above data sets is done with both 
DEA and SDA techniques. These methods are based on 
the prescription that the numbers in a time series {^i} 
are the fluctuations of a diffusion trajectory; see Refs. 
0,0,0] for details. Therefore, we shift our attention from 
the time series {^i} to the probability density function 
(pdf) p(x, t) of the corresponding diffusion process. Here 
X denotes the variable collecting the fluctuations and is 
refered to as the diffusion variable. The scaling property 
of p{x, t) takes the form 



p{x,t) 



(1) 



where 5 is the scaling exponent. The DEA (3 is based 
on the evaluation of the Shannon entropy, S{t), using the 
pdf If the scaling condition of Eq. holds true, it 
is easy to prove that the entropy is 



S{t) 



p{x, t) \n[p[x, t)]^A + 6 \n{t) , (2) 



where A is a constant. The SDA 0, instead, is based on 
the evaluation of the standard deviation D{t) using the 
same pdf (Q) and yields to 



D{t) = \J {x'^]t) - {x;tf oc t" , 



(3) 



where H is the Hurst exponent 0, E|- 

The Levy- walk statistics is a particular form of anoma- 
lous diffusion obtained by generalizing the Central Limit 
Theorem JTtI- While all theoretical details can be found 
in Refs. jjHS 1^ and in the enclosed references, here we 
show the properties of a Levy-walk using a model that 
simulates solar flare intermittency. We generate an ar- 
tificial sequence {tj } distributed according to an iplpdf, 



iP{t) (X 



1 



(T + 



(4) 



where T is a positive constant. In the computer- 
generated example, we use T = 90 and /j, = 2.5, that 
is, we use 2 < /x < 3, so as to violate the constraint of a 
finite variance on the central limit theorem that is nec- 
essary for obtaining Levy-walk statistics. The sequence 
of waiting times {tj} is used to generate a time series 
{^i} given by = 1 only at the occurrence of the events, 
that is, for i — [J2jLi'''j]i M the integer part of a. 
However, in the interval between two events, that is, for 
[SjLi '^il^ * < [J2Y=i 'Ol' assume (,i = 0. As shown 
in Ref. |1, 0, i] by using results of Refs. [H El [13 
the time series {^i} generates a diffusion process with a 
scaling exponent S = l/(/i — 1). Moreover, this particular 
approach to Levy statistics yields a diffusion process with 
finite second moments, as in the Levy- walk Conse- 
quently, it is possible to evaluate the scaling exponent of 
the standard deviation, H, via Eq. Fig. 2a shows 

the first 1650 data points in the computer-generated se- 
quence which correspond to the first 5.5 time units, 
of a 510,000 long data sequence {^i}. In Fig. 2b we show 
the sequence of frequencies {ft} of impulses for each time 
unit where we use 1 time unit — 300 natural units of t. 

The analysis is based on the comparison of the scaling 
exponents S and H measured by DEA and SDA respec- 
tively. It has been shown 0, that if a time series is 
characterized by Gaussian statistics, the two scaling ex- 
ponents are identical, that is, S = H. If, instead, the 
process under study is characterized by a Levy-walk, the 
two scaling exponents S and H can be both related to 
the exponent fj. of the waiting time iplpdf of the under- 
ling intermittent process generating the Levy statistics. 
In the Levy-walk case the scaling exponent H is slightly 
larger than S and the two exponents are related to fi and 
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FIG. 3: DEA and SDA of the computer-generated Levy-walk 
sequences shown in Figs. 2a and 2b. The two straight lines 
show the two theoretical scaling coefficients 5 = 0.67 and 
H = 0.75 for n = 2.5. 
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FIG. 4: DEA and SDA of the North-South hemisphere and 
Land- Ocean region temperature anomalies. The scaling ex- 
ponent H and 5 are reported in the figure, the error is ±0.02. 



one another via the following relation 0, 



1 



3-2H /i - 1 



(5) 



We refer to Q as the Levy-walk diffusion relation. The 
assessment of the property expressed by Eq. ((SJ, that 
is, the comparison between the value of the scaling ex- 
ponents S and H measured for the earth's temperature 
data sets and their relation with the scaling exponents 
d and fi measured for the hard x-ray solar flares waiting 
time series 4], is used in this letter to verify the con- 
jecture that the earth's global temperature inherits the 
intermittent Levy-like nature of solar flares, within the 
temporal range from few weeks to few months. 

In the computer-generated model we interpret the fluc- 
tuations in the frequency of the pulses shown in Fig. 2b 
as inducing a similar memory pattern in the monthly 
mean global temperature fluctuations. Infact, a higher 
solar flare frequency implies a higher monthly solar activ- 
ity. The DEA and SDA of the computer-generated data 
are shown in Fig. 3. These data are fit with Eqs. 
and (j3J yielding the two scaling coefficients S = 0.67 and 
H = 0.75. These coefficients are consistent with fj, — 2.5 
obtained by using Eq. (|3J); the known input value for 
this realization of the model. The large picture in Fig. 3 
refers to the analysis of the original 510,000 data-point 
sequence {^i}- The inset refers to the analysis of 1700 
data-point frequency sequence {ft} shown in Fig. 2b, a 
length compatible with the length of the mean monthly 
temperature datasets here analyzed. The two graphs in 
Fig. 3 show very good agreement with Levy-walk theory, 
that is, with Eq. because both time series produce 
similar diffusion trajectories. In particular, the small in- 
set in Fig. 3 shows the typical bifurcation associated to 
the two scaling laws that characterize Levy-walk statis- 
tics and that the statistics generated by only N = 1700 
frequency data points is rich enough to obtain a satisfac- 



tory pdf with Levy-walk scaling properties. 

Fig. 4 shows the numerical results by applying DEA 
and SDA to the temperature fluctuation data of the 
four regions of the earth. On the ordinate axis we plot 
\n[D{t)/D{l)] and 3(1)- S{1); thus, aU curves start from 
0. The scaling exponents S and H are obtained by fitting 
the first 10 points. This analysis for the Global tempera- 
ture (not shown in the figure) yields S = 0.90 ± 0.02 and 
H = 0.95 ± 0.02. The high values of the exponents im- 
ply a strong persistence in the temperature fluctuations. 
Both scaling exponents S and H are slightly larger for the 
South hemisphere than for the North hemisphere and for 
the Ocean than for the Land. A reasonable explanation 
is that the Ocean is an almost stationary system with a 
very high effective heat capacity. The Land, instead, is 
subject to stronger random fluctuations (as Fig. 1 shows) 
due to a lower effective heat capacity and a higher mor- 
phological variability due, for example, to the presence 
of deserts, forests, mountains and valleys. Finally, we 
note that the standard deviation scaling exponents H 
are larger than the diffusion entropy scaling exponents 
S and seem to fulflU the Levy-walk diffusion relation (O 
within the accuracy of our statistical analysis. 

The scaling behavior in the temperature anomalies can 
be related to the intermittency of the solar flare activity 
through the mechanism suggested by the model depicted 
in Figs. 2 and 3. In fact, the waiting time distribution 
of solar flares is characterized by an iplpdf of the type 
of Eq. (gl) with /i = 2.14 ± 0.05 0. This value of 
would imply a Levy-walk with S = 0.88 ± 0.02 and H — 
0.93 ± 0.02 (O for a shuffled dataset that would destroy 
any temporal correlation among flares and conserve only 
the Levy component; see Ref. about the effect of the 
shuffling on these data. We observe that the difference 
of the above two scaling exponent values gives H — S = 
0.05 ± 0.02, that is a value compatible with the difference 
between H and S measured for all the temperature data 
sets. In particular, H -6 = 0.05 ± 0.02 for the Global 
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FIG. 5: DEA of the unshufHed waiting time of solar flares se- 
quence during the period 1991-2000. The value of the scaling 
exponent is 5 = 0.94 ± 0.02. 



temperature and for both South and Ocean regions and 
H - 6 = 0.04 ± 0.02 for North and Land regions. 

To understand the actual values that we measure for H 
and S we first notice that the waiting time sequence of so- 
lar flares is weakly correlated 4J. Therefore, the effective 
value of 6 for the original unshuffled sequence is slightly 
larger than that calculated via the waiting time iplpdf ex- 
ponent fi for flares (5 = 0.94 ± 0.02, as Fig. 5 shows. 
We observe that this value of S is perfectly compatible 
with 6 — 0.93 ± 0.02 measured for the Ocean tempera- 
ture. The Ocean, in fact, is the region of the earth that 
should best mirror solar activity because it is the largest 



and most homogeneous region of the earth with a very 
high effective heat capacity; properties that retain corre- 
lated events in the temperature. Also we observe that the 
occurrence of uncorrelated events would slightly reduce 
the Levy-like memory and increase the Gaussianity of the 
data. This might explain why the scaling exponents H 
and 5 as well as their difference slightly decreases (0.04 
against 0.05) between the North and Land regions. Fi- 
nally, Fig. 5 shows that the scaling properties for the 
waiting time of solar flares last at least 200 consecutive 
flares and the waiting time for clusters of 200 solar flares 
may last for several months. The largest waiting time 
between two consecutive solar flares, during the period 
1991-2000, is almost 3 months and on average there 
are almost 60 solar flares per month. Therefore, the tem- 
poral resolution of these clusters is compatible with the 
temporal range of scaling of the temperature data sets 
that last at least 1 to almost 20 months, as Fig. 4 shows. 

In conclusion, the affinity of the scaling exponents ob- 
tained through our analysis suggest that the earth's tem- 
perature anomalies inherits a Levy-walk memory compo- 
nent from the intermittency of solar flares. We obtain 
reliable scaling properties within a short time interval, 
ranging from a few weeks to a few months. Morever, the 
joint use of SDA and DEA has proved to be very useful 
in detecting a significant Levy component in the time se- 
ries and, in general, it may suggest the presence of a link 
between a phenomenon under study and an intermittent 
process characterized by a waiting time iplpdf. 
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